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abstract 


Supervised Spectral Pattern Recognition and Visual 
Interpretation of MSS Imagery are two parts that constitute this 
thesis work. In pattern recognition, soils and surface features 
of the Krishna -Goda vary delta bearing orbit index No, l42--»049 are 
classified using K ^Class algorithm. Various combinations of the 
4 MSS bands are tried to achieve the best results » It is observed 
that a combination of more than two MSS bands gave very unsatisfac 
tory results. Further, a line printer rnap based oh K .Class 
Classifier decisions is prepared showing three- surface features 
viz, agricultural, waterbodies and built-up area, Ihe area 
occupied by each of the above classes is precisely determined. 

In the latter part, keys are developed for visual inter ■ 
rotation of MSS imagery to prepare soil-classif ication map/ built 
up area map and surface drainage map. It also explores the possib 
ilitios for salinity survey, turbidity studies and study of river 
a--,; irics using satellite imageries. Also, the lineaments preseni 
in the area and a Rose diagram to find their frequency azimuthal 
variations are drawn. For mapping purposes using visual inter - 
pretation techniques. Banc! .2 imagery (0.5 • o.6/w m) is found to be 

extremely useful. 
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CHAPTER X 

introduction AND OBJECTIVES 


1 »1 General 

October 1957 marked the beginning of the space age with 
the launch of SPUTNIK by USSR. For the first time in history, a 
man-made object was circling the earth, covering different regions 
and countries of the world from its vantage point in space. The 
years following the launch of Sputnik were the Cold War years 
with each super power tyying to dominate the other in space. This 
led to rapid development of satellite borne sensor for reconnains** 
ance purposes. Today, one has lost count of satellites that spy 
the earth every moment for military and civilian purposes* 

The civilian remote sensing programme was essentially an 
US and USSR effort in the seventies. U.S.A. made the data from 
their satellites available to many countries around the world at 
extremely reasonable costs. As a consequence the use of such 
data has become common the world over. With the advent of high 
resolution satellites like French SPOT System, remote sensing 
techniques are being reliably applied in many a number of fields* 

A number of studies like crop pattern determination, crop 
estimation/ landuse, surface water distribution, river course 
monitoring, forestry planning, geological mapping, fishing habitat 
characterstics etc. are being taken up* 
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Here in India, we have various facilities like data 
products laboratory at Space Applications Centre (SAC), full- 
fledged acquisition and processing facility at NRSA, Hyderabad for 
undertaking various projects for diverse applications. The launch 
of Indian Remote Sensing Satellite (IRS-1), no doubt, gives the 
indigenous efforts a further boost in this field. 

1.2 Satellite — its products 

The data obtained from satellites is converted into a 
variety of data products such as high density digital tapes (HBDi), 
70 mm film, microfische, black and white as well as-color prints, 
computer Compatible Tapes (CCT) and false color composites (ICC) 
by four different levels of processing. At first level browse 
products are generated in the form of HDDT and film negatives for 
all the bands after eliminating the cloud covered areas through 
quick look data. This product will be corrected for radiometric 
and earth rotation effects, annotated with salient details. The 
standard products are generated at level -2, that are corrected 
for sensor scene and platform-related geometric effects. Precision 
products are gere rated at level-3, and special products at level-4. 
Special products like floppy diskettes, use standard products 
(such as CCT) as inputs and are generated for specialized user 


s t 


needs for specific applications. 

The data products systems include image processing computer 
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special photographic laboratories equipped with s /stems for 
processing, developing, and printing of both B & W and color 
photographs and sophisticated recorders like laser beam recorder, 

1*3 Applications — - techniques 

Remote Sensing has proved its worth in a myriad of 
applications. Indian Scientists at various well equipped research 
organisations and frontline academic Institutes have applied the 
satellite products for diverse applications geological, agricul- 
tural, landuse , animal husbandry with diverse techniques— computer 
-aided,, visual and image processing , 

Snow melt run-off studies (Ramarnoorthi, 1983), ground 
water table prediction studies (Rampal, 1984) , soil mapping 
(Karale et.al, 1983), forest survey and management (Madhavan Unni, 
1983), land evaluation and classification for agriculture (Murthy, 
Venkataratnam and Saxena, 1983) are few to quote among Indian 
Works. Deekshatalu and Ktishanan, 1983 presents an overview the 
basic research problems in remote sensing in a discipline oriented 
approach. 

1 .4 Acquisition of Satellite Data 

MSS data of Landsat-4 has been acquired for the present 
study. Bandsat-4 was launched on July 16, 1982 and carried the 
Thematic Mapper (TM) and Multi Spectral Scanner (MSS) # The TM 
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collects radiometric data in seven spectral bands and has a ground 
pixel resolution of 30 m in the non-therrnal bands (120 m in the 
thermal band) compared to the four spectral bands of MS.S and low 
resolution of 80 m. Landsat *-4 TM data cannot be acquired at the 
Indian Landsat Earth Station (ILES), situated at about 60 1<M 
from Hyderabad/ A.P., due to the failure of X— band radio channel 
in Feb/ 1983 which is responsible to transmit the Til data directly 
to ground stations. The MSS r however,/ has not been affected as 
it uses a seperate S«band radio channel. Landsat—5/ launched on 
March 1/ 1984, supplies both MSS and TM data. This space craft is 
positioned on the World Reference System (WRS) with its cycle 
offset 8 days from that of Landsht «• 4. 

The attitude of Lands at 4 and 5 (705 I<M) compared to 
that of Landsat 1/ 2 ard 3 (9l0 KM) means that a different WRS 
path-row number is used to refer the nominal scene centre. While 
the row number remain same, path numbers are 10 to 12 lesser fehan 
those for the first three satellites. The MSS bands of land sat 1, 

2 and 3 are numbered 4, 5/ 6 and 7 while those of the latter two 
are 1, 2 / 3 and 4. The TM spectral bands are numbered from 1 to 7, 
The RBV subscenes are numbered as A, B, C and D. The spectral 
ranges of these bands, resolutions and their major applications 
are listed in Table 1.1. 

The products obtained for the study are of the scene with 
index number 142-049. a CCT, black and white paper prints, films 
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of Bands 1, 2 and 4 are obtained from NRSA* Band-3 print, and 
film were not obtained because of non-availability. 

The video data obtained on CCT is in Binary mode/ 32- bit 
format which is not Dec system— 1090 Compatible. Dec is a 36— bit 
Word machine . The obtained tape is formated for 32-bit word 
machines like Vax, IBM and PDP series. If the tape is tried on 
Dec— lo, it tries to read the four bits of the next word and thus 
tampering the data. As a remedy to this problem* a dummy blank 
was padded after 4th, 9th, 14th..,. bytes resulting in the 
increase of the record size from 3596 to 4500 bytes. This was 
carried out by Nagaraju (1986) with the helip of CMC people. This 
conversion can be carried out with DEC-10 also with the help of 
CHANGE- •EXE, a system program. This program is given in the 
Appendix-1 succeeded by various switches that may be used in the 
program. It should be borne in mincl that DEC-lo cannot support 
800 BPI tapes and so they cannot be tried for a change in their 
format, DEC— lo which accomodates only 1600 BPI tapes reads the 
800 BPI tapes at a speed twice of what is desired. NRSA supplies 
both varieties of tepes and 1600 BPI tapes only may be acquired 
for the local use. 
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1.5 Objectives of the Study 


The following objectives are envisaged in the present 
study , the study area being Krishna -Godavari delta having an 
orbit index No. 142 049, Two techniques are employed to 

achieve these objectives viz* supervised spectral pattern 
recognition (K-Class Classifier) and visual interpretation 
techniques using MSS data and paper prints. 


With K- -Class Classifier using the signals of four MSS bands 

(1) to classify the soils of the region 

(2) to classify the surface features such as water bodies, 
agricultural area, built-up area etc. 

(3) to prepare a line printer map based on the classifier 
decisions showing the surface features. 


( 1 ) 

( 2 ) 

(3) 

(4) 


With visual interpretation o£ three MSS prints: 

to prepare a soil classif ieHtion map 
to delimatu the built-up area in the delta zone 
to update the drainage map of the area 
to draw the lineament map of the study area. 
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CHAPTER II 

PATTERN RECOGNITION USING K-CLASS CLASSIFIER 

2*1 P attern Recognition—' Various Stages and Techniques 

Spectral Pattern Recognition is a computer aided technique 
by which MSS digital data may be analyzed quantitatively and 
automatically. The advantage with this analysis is that digital 
data of desired bands may be used in any combination. In this 
process , we use the computer to look at the multiple channels of 
digital data. By dealing with the image data quantitatively/ the 
spectral Information of any number of channels can be fully 
evaluated* 

The concept of representing MSS data in a numerical format 
is illustrated in Fig. 2.1. The figure shows MSS measurements of 
certain features in a landscape made along one scan line* The 
digital values or spectral responses of these features in five 
bands arc represented with vertical bars* If these spectral 
patterns are sufficiently unique for each feature typO/ they may 
form the basis for an automated interpretation of the image data 
using a spectral pattern recognition proccedure* There are 
various techniques available to analyse and classify a set of data 
into groups or classes* These techniques can be broadly divided 
as supervised and unsupervised pattern recognition methods* The 
difference between the two is that the former needs a training set 




■ ioJiUwsng spec 
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of data acquired with human knowledge. 

In unsupervised classification methods the mode seeking 
alogorithm is employed to classify the data into groups of points 
or clusters of similar spectral characteristics f Hence collection 
of representative training sets is circumvented. In class-room 
anology, these methods are learning with and without teacher. 

The three basic steps involved in typical supervised pattern 
recognition procedure are 1) training stage , 2) classif ication 
stage, 3) output stage. Ihis is illustrated in Figure 2,2. 

The training stage, a requisite process in supervised 
classification, involves compilation of an ‘interpretation kcy J for 
each feature type of analyst^ interest. 

To arrive at the key, the analyst should have the knowledge 
of the “training areas” or “training sites" wherein the features 
of one's interest are contained. These may bo termed, as reprcsentn~ 
tivc sample sites of known cover type. The digital data at these 
areas are retrieved from Computer Compatible Tapes (CCT's) to 
form the training set. 

In the classification stage, each pixel in the image data 
set is compared to each category in the numerical interpretation 
key. This comparison is made numerically, using any one of a number 
of different stratigies to decide which category an unknown pixel 


* 


value "looks most$ like" 
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Each pixel Is then labeled with the name of the category 
it resembles/ or labeled "unknown" if sufficiently similar to any 
category. The category label assigned to each pixel is then 
recorded in the corresponding cell in an interpreted data set-. 

Thus, multi-dimensional digital data is used to categorise the 
area of interest. After the entire data has been categorized, 
the results will be presented in the output stage, commonly in the 
form of a map. this way, we are also able to enlarge the imagery 
without the loss of detail or contianity. 

Because it is numerically based, spectral pattern recog-' 
nition is largely an automated process. In this respect, visual 
image interpretation and spectral pattern recognition are 
complimentary procedures* 

Another classification of pattern recognition techniques 
is parametric and n on-parametric methods. In parametric methods 

iU.A 

each pattern class is characterised by a statical distribution 
which in turn is dependent on certain number of parameters. The 
non-parametria methods do not assume any such distribution. Tabic 

2.1 gives a partial list of classification methods. Two of the 
methods are briefly described here. 

2.2 ( a ) Ba yes 1 c lasslfier 

One of the most widely used parametric classifiers is 

In parametric classification a probability 


the Bayes 1 Classifier. 
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density function is assumed and the parameters' of that distribution 

aie estimated. let ,. #/ be random variables where 

x i t ^ riG nois y Measurement® of the ith feature. Let p(x/C.) be 

J 

the conditional probability density function of class j and p(Cj) 
is the a priori. probability of class Cj, The task of the classifier 
is to assign the input sample such that the probability of misrec*- 
ognition is minimized. 

The Bayes 1 decision is that 


X un C 


if 


P (C i )p (X/G i ); > ^•(Cj)p(X/Cj) / for all j 


Assuming Gaussian distribution with moan vector and Covariance 

matrix k; 

i 


P (X/CJ « 


1 ( 2 A 


exp 


-Jj (X-M i ) T K i “ 1 (X - M, j 


then the decision boundary between classes i and j becomes 


log 


p(Cj_) 

Flop h 




T., -1 


(XmM ± ) K ± ( x *^ i ) - (X^j) Kj (X-Mj) 


= 0 


The above rule is also referred to as the maximum likelihood 
classification estimation rule (MLE) which has been very popular 
for classifying remotely sensed data. 
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2*2 (b) K-Class Classifier 

The K^Class Classification algorithm developed by 
Zagalsky (1968) is a supervised and non-par ametric pattern recog*- 
nit ion method. It is based on the derivation of a transformation 
matrix (B). The main ccmputing steps are given below. For 
mathematics and algorithm of the classifier, Serreyn and Nelson 
(1973), and Sashi Kumar (1982) may be referred to; 


Computing Steps : 


a ) 


Compute the a priori probability p 
class i where 


„ No of points in class i a . 
Total no. of points 


. of occurrence of each 


— 1 ^ 2, 3, #»«... K< j C las se s 


b) Calculate the mean of the attributes (signals or grey levels 

retrieved from CCT) over all classes. 

K<"(X. of class i) (No, of points in mode i) 

5? . = _ 4 — 1- — 


1=1 


Total no, of points 


for j - 1, 2/ ***.*.../ N features 


G ) 

d) 

e ) 


Calculate the attributes covariance matrix 
T T 

0 = L x Y - x x J 


*•1 

Calculate the inverse of Covariance matrix (|2P ) 
Calculate the Transformation matrix (B) 
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Calculate the Vector of Constants (C) 



g) Calculate the elements cL of the class vector d for the 

attribute vector X , where 

d i 88 L BX * c i + 1 J P± ? i - 1/ 2, * . . , . K class 

h ) Assign the attribute X^ to class i for which d^ is maximum 
This classifier as a computer program is implemented to 
various applications in this study. To train the K*-Class Classifier 
the investigator needs to have a training set of data., This set 
consists of sample data from each class of interest. Then all 
the above steps are performed, The only unknown in the equations 
for the decision vector elements is thefeature vector to be 
classified. The decision is calculated by selecting the maximum 
element value dj. The feature vector X is assigned to Class j. 

By this procedure, the map of an area can be made representing the 
classes by different distinct characters. 

A brief overview of pattern recognition and image processing 
is presented by Deekshatulu and Kamat (1984) and Pu (1984), One 
may refer to these papers for extensive coverage of pattern 
recognition techniques. 
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2 * 3 gggv lous Work With the Classifier 

Serreyn and Nelson (1973) applied this technique on data 
taken from ERTS imagery for classification of follow, Corn and 
Soyabeans* Lidster et al,(l9 ) applied multiple regression, mot]e 
seeking, and K-Class classification analyses for correlating 
digital data of bandsat and Aircraft imagery with Water table 
depths in irrigated agriculture. The mode seeking, and K*-Class 
classification analysis of a Gorn field resulted in the correct 
classification of 91 % of water tables. Sashi Kumar (1982), 
employed the algorithm to find the number of classes that occur in 
a two dimensional classification with the groundwater table depth 
and the Band - 7 spectral reflectance for the alluvial portions of 
U.P. and Rajasthan , 

2 • 4 Applications of the K^Class Classifier 

The classifier is put to variety of uses in the present 
work such as soil identification, surface feature identification, 
and mapping with the help of classifier. Each of these applica- 
tions is dealt in detail in subsequent sub-sections. Various 
c combinations of bands are tested to see how they respond to the 
classifier, 

2.4.1 Soil Identification 

The study area consists mainly throe- typos of soils viz;-Recen1 




alluvium (sedimentary unconsolidated) sandstone (sedimentary 
consolidated ) and unclassified crystalline rocks as it turned out 
from visual interpretation of satellite imagery. The classes 
chosen for this ore Recent alluvium (black soil of delta area), 
unclassified crystalline rocks , Khondalitcs, Coastal Sandy Soils, 
and Char nockites , The exact locations of the occurrance of these 
soils are noted from survey of India maps. The geographical 
Co-ordinates i.e. latitudes and longitudes of these points arc 
converted into scan line numbers and pixel numbers. 

The basis for the selection of the above variables ol 
soils is Khondalite and Charnockite both being crystalline fall in 
the zone of unclassified crystalline rocks* These two classes 
are to be shown as seperate entities after classification if and 
only if the place where attributes were obtained for unclassif iec. 
crystalline rock class contains a different material i.c. other 
than the above two. 

To start with , a 3-class, 2 -feature problem is chosen* 

The term feature here means the data from a band of MSS. The Bands 
selected are 1 and 2 of MSS (Landsat 4) and the classes are 
Khandalites, unclassifieds, and Recent alluvium. In each class, 
eight attributes (pixel values) are rotrived and then subjected, 
to K-Ciass classification. The results are produced as Table 2.2 
The number of signals being same in each class the probability of 
occurrance of each is given as 0.333* The transf ormation matrix 
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is represented by B of size (3x2). From the table of decision 
vectors it may be observed that the value for whichever class is 
maximum, the decision will be that the signal is attributed to 
that class. Thus r for the sample No.l, the decision values are 
1.09, o,07 , and -0.16. So, the sample j.s placed in the class 1 
thus making the- correct identification. The confusion matrix shows 
that no, of miscalculations is zero (the diagonal Elements are 8 
each in value ) . The decision boundaries are shown next in a 
two-dimensional space. 

Next, the attributes of the same three classes are retrived 
from Band-2 and Band-3 of MSS to put before the classifier. This 
time also, the classifier has successfully identified the samples 
with lOO % accuracy. The results are in Table 2.3. 

Having satisfied with the results with the help of 2 —band 
signals, it is tried to test the sarre classes with the help of 
Bands 2, 3’ and 4 signals. Thus it becomes a 3-class, 3-Feature 
problem. In each class, eight signals were given as samples to 
the classifier. Now, it is interesting to see that percent 
correctly identified has fallen down to 66.67 %. Earlier, in 
both cases, it was 100 The confusion matrix (Table 2.4) shows 
that the four signals of class (3) are identified with class (1) 
and the next half with class (2) resulting in a drastic fall of 
accuring. The examination of the spectral responses of these 
classes in the three bands may give an explanation for this fall. 
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Be t us consider the response of Khandalites in Band-2, These 
range from 21—25 very much similar to that of Recent alluvium 
(20-26) in same band. Also, the grey levels of unclassifieds in 
Band - 4 (45-61 ) is also similar to those of alluvium of same 
band (59-71), Coming to the response of alluvium in Band 3, some 
of them are similar to those of class (2) of same band. Thus, 
the third class lost its ’identity’ during the calssification stage* 
An inference drawn from this is that the bands we choose for various 
classes should give unique identity to them. This may be achieved 
with a data of good number of bands, However, it is found later 
that a combination of two bands gives excellent results for 
mapping with the classifier. 

Next, another class is added to these three and tested as 
4-Class, 2-Peature problems. The classes are Khandalites, 
unclassifieds. Recent alluvium and coastal sandy soil and the 
features are Band —1 and Band-3 responses. It is worth noting 
that the classifier give 100 % accurate results in this case also. 
With the same bands, charnockite as the fifth class, the classifier 
identified the classes with 72 % accuracy. This may be due to the 
association of charnockites with unclassifieds. The outputs of 
these two classifications and co-ordinates of the classes on the 
imagery are not available for reproduction because of system H/I 
problems at the last phase of the work* 

next step, charnockites has been dropped and a 


So, in the 
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surface feature (water bodies) is substituted. The results shoved 
that two of the unclassifieds are identified with sandy soils and 
the rest with water bodies. Thus / the whole class is misidentif led 
by the c las sifier (Table 2.5 ) % 

2.4.2 Surface Feature Identification 

The study area includes many interesting surface features 
such as rivers, reserve forests, lakes, marshy lands, coconut 
plantitions along the sea coast, built up areas. Location of 
these features is made with the help of topographical maps of the 
area. These geographical co-ordinates are converted into image 
co-ordinates to enable the data retrieval from CCT. Three of these 
features namely water bodies, built up areas and dense jungles are 
chosen for the computerised statetical classification. Dense 
jungles are Muttayyapalem Reserve Forests (RF), Water bodies 
(Krishna river and Bay of Bengal}, built up areas being urban regions 
of the area. Table 2.6 gives the geographical and image coordinates 
of these locations. 

First, eight samples each of water bodies and built up areas 
ore subjected to classification and the outcome is free of any 
miscalculations (Table 2.7). Bands used for this are 3 and 4 of 
MSS. Next, another 2-class, 2-features problem is developed with 
water bodies and dense jungles as classes and Band-3 and Band. 4 
values as features. This time also, the results are 100% 
accurate (Table 2.8). 
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7 -CUSS# 2-Ff-ATURE PROBLEM 


t K classes are i ' * 

nonius and built-up areas in band-3 & sand-4 


N " L A S 5 a 2 r/FfcJArx 2 NSTGb 16 NPRl)B« 1 


ru?: 

S l 

SNALS IN C 

LASS 1 APE 

I 




2B.0J 

1 1 .00 





31.00 

B ,00 





it .oo 

7,00 





34.00 

11.00 





20.0 0 

11,00 





24.00 

7.00 





22.01) 

11 .00 





36.00 

n .oo 


VHP 

si 

GNALS IN C 

LASS 2 ARE 

l 




as. og 
97.00 

53.00 

61.00 





77.00 

63.00 





77.00 

63.00 





91,00 

67.00 





70.00 

57.00 





9R.00 

77,00 





96,00 

65.00 


number 

nr 

SIGNALS 

IN CLASS 

Cl) 

Numbs 

R 

OF 

SIGNALS 

IN CLASS 

(2) 


TUE matrix b IS I 

-0.01225 ^O'aO‘2 41 1 

.0.01275 0.02411 

DECISION VECTORS 


SAMPLE, no; class 



DU) 

D (2 ) 

0.97 

0,03 

0.97 

0,03 

1,00 

-0.00 

0.94 

0.06 

1 .02 

-0.02 

t.05 

-0,05 

1.01 

-0,01 

0,93 

0,07 

0.11 

0,89 

0,04 

ojg 

0,05 

0.95 

0,05 

0.95 

0,09 

1.09 

0.16 

0,84 

0,19 

1,19 

0,03 

1,03 


J.— j. W„ 

DECISION 



2 


J'-' 


CLASSIFIED ?\S. 


CI.AflS 1 2 

1 fl 0 

2 0 8 

PRS&A 6 [LT TIES* 0,500 0 . 5 OQ 

Mf). j r Mis:AiiCui,noMs *= o 


pr.rcfiNT roff«r:rLr EUG.mriF:o over alls t oo.ooooo 

pEn:?NTAlJE MATRIX FOR ABOVE CLASSIFICATION 
CLASS 1 2 

1 100,00 0.00 

2 0.00 100.00 

> 

DECISION BOUNDARIES 


the bound A ny BETWEEN CLASS 1 AMD CLASS 2 IS ; 

-0.01225063 XI + -0.024U196 X2 a -1.556959 20 

THE BOUNDARY BETWEEN CLASS 2 AMO CLASS 1 IS I 

0.01225063 XI + 0,02411196 X2 s 1.55695820 


TABI£ 2,8 
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2-CMSSr 2»FP:>\TUnK PRJALEH 
V'.\r. CLASSES ARE t 

4 A TER ROMES AND DFNsL JUNGLES OF* BANO-3 $t RAND-'l 
CLASS* l N FEAT* 2 NSTGs 16 MPPQB* 1 


r h 

SIGNALS IN 

CLASS 1 ARK 1 


28.01) 

It .00 


l* .00 

ft , 0 0 


3 t .00 

7.00 


31.00 

It .00 * 


20.00 

11 .00 


24.00 

7.00 


22,00 

It .00 


36,00 

11.00 

THE 

SIGNALS in 

CLASS 2 ARK l 


81.00 

62.00 


93,00 

B7.00 


110,00 

98.00 


77.00 

62.00 


7-1.00 

59.00 


37.00 

71*00 


35.00 

70. 00 


81,00 

63.00 


N HMB £ R 3F SIGNALS IN CLASS (!*)» 8 
JJdMBEfl OF 'SIGNALS IN CLASS <2)« 8 

rffE 'HT.RI)C l n is s 

-0,01726 -0.01387 

0.01726 0.01387 

'■ DECISION vectors 

rm „ M c M M MMW.W, •(»'*» 


SMPLE 00. class Ml) D ( 2) DECISION 




CLASSIFIED AS. 


*•*»*»/• 

CLASS 1 a 

1 9 0 

2 0 8 

PROBABILITIES* 0,500 0,500 

NO.’ OF* MISCALCUL TIQNS » 0 


percent correctly identified overalls* ioo.oooqo 


PERCENTAGE MATRIX FOR ABOVE CLASSIFICATION 
CLASS l 2 

1 100,00 0,00 

2 0.00 100,00 


DECISION BOUNDARIES 


THE BOUNDARY. BETWEEN CLASS 1 AND CLASS 2 IS j 

-Of 01725602 XI + -0,01 387335 X2 * -1,55172370 

THE BOUNDARY BETWEEN CLASS 2 AMD CLASS 1 IS t 

0.01725602 XI + 0,01387335 X2 * 1,55172370 
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Having known that the classifier is able to classify two 
classes at a time successfully, now, all the three (water bodies, 
built up areas and dense jungles in order) of bands 3 and 4. arc tried* 
After classif ication, it is found that four of the built up area 
samples are identified with dense jungles while two of the latter 
class are identified with built up area class. The percent 
correctly identified overall is 75.0 % (Table 2.9). To improve 
this percentage, the combination of Band— 2 and Band -3 signals is 
tried (Table 2„l0l. This time, all the eight signals of the third 
class are correctly identified while three of built up class are 
recognised with third class. The percentage of correct identifi- 
cation, thus, bee ones 87.5 %,a marked 12.5 % increase over the 
previous combination. Other combinations may also be tried to 
see whether it still can be improved. 

Finally, a 3 class- 3 feature problem Is developed with the 
some classes of bands 2, 3 and 4 signals (Table SL.ll ) „ The 
percentage of correct identification, expectedly., fell down to 
54* 16 % . The confusion matrix shows that three of the first 
class samples associated with the second one and the second class 
samples totally misrepresented. 

2.4.3 Mapping with the Classifier 

Now that we are sure that the classifier gives more than 
satisfactory results with a combination of 2 bands, a small region 


TABLE 2.9 


3 -CLASS, 2 -FEATURE PROBLEM 


the 21* a sses a he t . . „ „ 

WAVER BODIES, BUILT-UP AREAS AMD DENSE JUNGLES OF BAND-3 t 


bamu-4- 


M-LASSb j NFEAr* 2 WSIG« 24 MPRQfla 
THE SIGNALS t N CLASS 1 ARE I 


28.00 

34.00 

31.00 

34.00 

20.0 0 

24.00 

22, 0J 

36.00 


11.00 
0.00 
7.00 
1 1.00 
it. oo 

7,00 
11 ,00 
11*00 


the signals in Class 2 are 1 


86.00 

82.00 
77.00 

77.00 

91.00 

70.00 

08.00 

86,00 


53.00 
61*00 

63.00 
63*00 

67.00 

57.00 

77.00 

65.00 


THE SIGNALS IN CLASS 3 ARE 1 


81,00 

93.00 

110.00 

77.00 

74.00 

87.00 

85.00 
81 .00 


62.00 

87.00 

98.00 

62.00 

59.00 

7 i.no 

70.00 

63.00 


NfJ.MBSR OP SIGNALS' IN CLASS (1) = 
N'TMRSR OP SIGNALS IN CLASS (2)* 

number op signals in class (3)* 
i'HE MATRIX fl is I 


- 0.03354 

3,07374 

- 0,04020 


- 0.01546 

- 0*05331 

0.06577 
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hfcXISIOU V FC TORS 


l M h*.» w «*|*,**^ M A* 0 *** •» W m ■! M ' 


SrtMPLS : 4 H . 


9 

10 

} 

H 

H 

15 

1 § 

17 

tfl 

19 


CUSS 


n( l ) 


0.91 

9.99 

0.93 

0.08 

1.03 

1.01 

1.01 

0.86 

0,09 

0,08 

0.13 

0.13 

■0.05 

0.24 

-0.07 

0,02 

0.09 

•0.17 

■0.42 


0(21 

DID 

0.0 \ 

0.02 

0,23 

-0.12 

0.17 

-0.10 

0,18 

-0.06 

0.16 

0.00 

-8:1 1 

0.12 

0.10 

0.23 

-0.08 

0.75 

0.17 

0.52 

0.40 

0.36 

0.51 

0.36 

0.51 

0,64 

0.41 

0.29 

0.47 

0.40 

0.67 

0,55 

0.43 

0,48 

0.43 

0.36 

0.82 

0.59 

0.83' 


DECISION 


ll 

l\ 

24 


1 0.1 3 0.3 8 0.4 8 

3 0.18 0.36 0.46 

3 -0,02 0.4R 0.53 

3 0.00 0.44 9.55 

3 0.08 0,46 0.45 

classified as 


CUSS 1 2 

1 8 0 

2 0 4 

3 , 0 (i> 

PROftA^rLUrKSa 0.3 33 0.3 33 

NO , 3 F MlSCAliCUoriCJMS a 6 . 

correctly identified overall* 75.00000 


3 

0 

(i 

6 

0 
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PSRCF'I f/\(?£ MATRIX FOR AhbVFT CLASS CriCATIOiN 
CLASS l 2 3 

1 100,00 0.00 0,00 

l 0,00 50.00 bn, 00 

J 0,00 25 ; un 76.00 



utCtmiih rum - j r* a r» i r a 


rn~ suomdary betwceh class j and class 2 is i 

-0.03575919 XI t 0,01161466 X7 s -1.7R625150 

Trip: B ~J U >1 0 A R Y 0EIWEEM CLASS 2 AMO CLASS 3 IS J 

t. 07*526? XI + -C ,C SBfcSOfi-b A5 a 0.629002 8 2 

. )Hf aot^o/py atTwepM ci .# ps z no cr#fs i ir ; 

-0.00222289 XI + 0. 02707600 X? * 1,15724070 


TABLE 2,10 
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3 'CLASS, 2-Kf?AT PROBLEM 

we glasses are i 

^ APKR BODIES, BUILT-UP AREAS AND DENSE JUNGLES OP B AND**2 h BANt)*3 

N:LA55S 3 NFEM’= 2 NST G= 24 NPROB= 1 
THE SIGNALS IN CLASS 1 ARE I 



52.00 

28.00 


56.00 

34.00 


52.00 

31 ,00 


58,00 

34.00 


39.00 

20,00 


28,00 

24.00 


29,00 

22,00 


50.00 

36.00 

the signals in :: 

LASS 2 ARE t 


71.00 

B6 .00 


37,00 

0 2.00 


68.00 

77.00 


24.00 

77,00 


50.00 

91 .00 


28,00 

7 0,00 


47,00 

86,00 


4 6,00 

92.00 

the signals in : 

LASS 3 ARE l 


39,00 

ouoq 


43.00 

93,00 


41,00 

110,00 


27.00 

77.00 


3 3,00 

74,00 


4 4,00 

07.00 


33.00 

85,00 


■ 34,00 

81 .00 

number 

3F SIGNALS 

■ IN CLASS (i)« 

number 

JF SIGNALS 

IN CLASS (2)3 

number 

OF SIGNALS 

In class (3)8 

f HE matrix n IS t 


.0 . 

01371 -0. 

04907 , 


02599 0. 

02237 

*0 , 

0 3 97 3 0. 

02670 



DECISION 


ijEt: rsi’i 


3MPI.3 N '.} • U.ASvS 


Dtl) 


0(2) 


on) 


5 

6 
7 
R 
9 

lt> 

U 

If 

1.7 
18 
1 9 


21 

24 


1 

1 

2 

2 

2 

2 

2 

2 

2 

2 

3t 

5 

3 


O.g? 
0.01 
0.94 
0,02 
1 .06 
0.05 
0 , 0 ti 
0.B5 
0,13 

o.o 4 
0.26 
0,06 
- 0.01 
0.10 
0,02 
-0,08 
O.Qb 
-o.ii 
- 0.00 


0.08 
0.1 5 

-o.oi 

-0,03 

0.04 


0 1 2 t 

o’.iS 

- 0 . 0 } 

-0.11 

- 0.12 

o.n 

0.73 

0.40 

o.&3 

0.25 

0,65 

0*23 

0.52 

o.:s5 

q .41 

0.53 

0.64 


0.28 

0.31 

0.50 

0.39 

0,37 


3,1 2 
1 D . 1 2 
’3.0) 
■0.14 
•0.02 
0.15 


0 

1:11 

m 

m 

5 2 

j',.5 9- 
0 .‘7 5 


1 


0.6 5 
0.54 
3.51 
0.64 
0.59 


I 


CLASSIFIED AS 


CLASS 1 2 3 

1 BOO 

2 0 5 3 

3 0 0 0' 

PR1BA'31IjI TIES- 0,3 33 0,333 0<*,3 3 3 

?I0,. 3? IISCALCULTICJNS - 3, 

CORRECTLY tOilMTlFl-ED OVERALL* 07,50000 

i 

PSRCSMFAGE MATRIX FOR ABnVE CLASSIFICATION 
CLASS 1 2 3 

t t 00,00 0.00 0,00 

2 0,00 62,50 37.50 

3 0,00 0,00 100,00 




T-HE 3. ■JSOARY HKI'rfSSW CLASS J AMD CLASS 

-a.oumis? xi ♦ -o.ommi x 2 . 

rHB 3: ”' SM *' f “SfWW CLASS 2 AMD CLASS 2 

0.0^18^379 xi t "0.U0144339 X2 a 

rue addmoarx between class 3 amd class 1 
* 3.01780222 XI f 0.02525530 X 2 a 


IS I 

-1 - U24226D 
IS I 

0. B51&9D85 
IS l 

0,99133175 


TABLE 2*11 


3»CLAS & t 3-KEATUPS l 

r hr classes ariT! 

ROOIES, BUILT-UP AREAS, and [)RNS£. JUNGLES IN 8-2, B-3, £ 8*4 


CL\SS= 3 NHCATb 3 NSTGs 

24 NPROB 

ritr SIGNALS T.N 

CLASS 1 ARE J 


52.00 

28,00 

11,00 

56.00 

3 4,00 

8,00 ’ 

52.00 

3l*»n) 

7.00 

58.00 

3 4,u 0 

11 ,00 

3 9.00 

2-0.00 

11.00 

28.00 

24, OD 

7,00 

29. DO 

22,00 

11,00 

5u • 00 

3 6,00 

11.00 

l'N[C SIGNALS IN 

class 2 are : 


71,00 

86.00 

53,00 

37,00 

82,00 

61 ,00 

68,00 

77.00 

58,00 

24,30 

7 7.00 

63,00 

5 8.00 

91.00 

67.00 

28.00 

70.00 

57,00 

47,00 

86,00 

65.00 

46,00 

92,00 

70,00 

rrtfi SIGNALS IN 

CLASS 3 ARE I 


39,00 

81,00 

62.00 

43,00 

93,00 

B7.00 

. 41,00 

lll),QO 

98,00 

27,00 

77,00 

62,00 

33,00 

7 4.00 

59,00 

4-1,00 

87,00 

71.00 

33,00 

05.00 

7U.OO 

3 4 , 0 0 

81.00 

63.00 


IMMBS'R of signals IN CLASS (1)s 8 
NiiviBEA JF SIGNALS IN CLASS ( 2 ) c 8 
N.HBER OF SIGNALS IN CLASS (S)s 8 


THE MATRIX 0 IS t 

0,62 63 0 -0,09400 0,00720 
-3,27690 0,0.351 7 -0,08067 
-0,3 3 940 0,05883 0.07348 


DLCJ-KiltU VECTOk* 


5A«fPli£ 9 3 . CLASS 0(1) t> C 2 5 


1 

1 

3.2b 

0,06 

2 

1 

3.90 

-0.15 

1 

i 

1:11 

-8:11 

5 

1 

0. RO 

1.17 

5 

1 

-1,63 

2 ,34 

7 

l 

-1,35 

2,12 

8 

1 

2.60 

0,34 

9 

5.51 

-2,14 

10 

\ 

-1,44 

0.74 

12 

\ 

5.18 

-3,99 

"?:I8 

H 

\ 

-1:18 

-1 a 26 
- MJ 

15 

2 

0.53 

-0,25 

16 

2 

0.1 5 

-0,22 

17 

3 

-0,99 

0.51 

lB 

3 

-0,4 7 

-0,39 

19 

3 

-t , 39 

-0.30 

20 

3 

-3.37 

1,58 

21 


-2,03 

1.07 

22 


-0.11 

-0.12 

23 

3 

-2.35 

0.90 

24 

3 

-2.03 

0.95 


classified as 

ASS J 2 3 

1 5 3 0 

2 3 0 5 

■3 0 0 8 

PRDfJABrta 1‘TfiSs 0,333 0,333 0,333 

M;>, OF 'IIJJC/UjCUL no NS all, 

KRCP.NT COUKECfLV IDENTIFIED OVERALL- 54,1 6667 



- 2 . '33 
-2.-75 

:i:$ 

‘Ml 

0.-2 3- 

-1 .94 
-2.37 
1 .>7 0 

“ 3 : 4 ! 

1 .47 

r .'B 6 
2.69 



2.-7 9 
! 

1.23 

2.43 

2 . o a 



PS.RCSN 1' AG E MATRIX FOR ABOVE CLASSIFICATION 
CLASS 1 2 3 

1 62, 5U 37,50 0-00 

2 37,50 0.00 62,50 

3 d.OO 0.00 fOO.OO 
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of the study area is chosen to generate a computer map based on 
K~Class Classifier decisions. The coordinates of the area are 
shown in figure 2.3. 

A(10 3 3, 701 ) B ( 1033 , 820 ) 


64 x 120 
matrix 


0(1096, 701) 0(1096,820) 


Fig. : 2.3 


Fran the toposheet, it is observed that there are only two classes 
present in the area viz water and non-water area. The criterion 
for the selection of this region is that in this small region/ 
Kommamur Canal and Nallamara Vagu intersect each other. Also, 
a reservoir is present near the end O, Tnc digital data of Band-2 
and Band-3 of this area is retrived for the analysis. The 64 x 120 
data matrix is assumed to contain tv;o classes and the computer map 
is generated (Fig. 2.4). 

It is worth noting the date of satellite pass for the data 
obtained ( 28th September 1983). This period is a harvest period 
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and the crop is generally paddy* So, the classes actually 
identified are agriculture and water bodies. The number of pixels 
present: in each class is given below the map. 

Now, some more- lines of data proceeding this area arc 
added to this making the co-ordinates of A and B as A{l00l, 701), 
5(1001, 701). This is done because this added area contains a 
village namely Return. The map generated by the classifier 
beautifully shows this village (Fig. 2.5). Another feature 
uncovered in toposheets is also observed on the map. This being 
a regular (hallow, aguare-like ) feature this is to be taken as a 
roan-made featuee. This is identified as water body but it may be 
deducted from its dimensions as a foundation trench filled with 
water. Number of pixels present in each class and area occupied 
by it in percentage is given as its annotation. 

The software developed for this work is explained under the 
head SOFTWARE in Appendix -2* The listings appear in subsequent 
appendices. 
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2.4.4 Observation? and Conclusions 

The classifier with 2-class and 2-band problem gave 
almost loo % correct identification. With 3-band signals of 
MSS, the percentage accuracy drastically falls down. The 
classifier gave 100 % accurate results in identifying Khondalities 
Recent alluvium. Coastal Sandy and unclassified soils. Infact, 
it can be deducted what exactly anclassif ieds are with different 
trials with the classifier. 

Coming to the classification of surface features, 2-band, 

2 -.class cases gave satisfactory results whereas while classifying 
waterbodies, built .up areas and dense jungles with Bands 3 & 4 
signals, it dropped down to 75 % . Ibis is improved by 12.5 % 
with Band*. 2 and Band-3 signals. 

The line printer map showing the three classes viz. 
agriculture, water bodies and built-up area prepared on the 
basis of classifier decisions matches with the topo sheet 
details. It, infact, showed some more details. 
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CHAPTER « I IX 

VISUAL INTERPRETATION OF MSS IMAGERY 

3.1 Introduction 

Though quite young , the field of aerial photography 
developed a number of broad divisions, each in itself having a 
remarkable diversity. Of its many offshoots, aerial photo— inter— 
pr station has acquired prominence because of its effectiveness of 
the technique. The technique includes the character sties like 
identif ication of man-made features, common terrain features, 
feature analysis^use of stereoscope etc. Though being an effective 
method, it is costly in view of the flights, man-hours and 
sophisticated equipment that be used in the process to get the 
photographs. Now that,' many satellites are operational and their 
products available at reasonable prices to the users round the 
world, visual interpretation of satellite imagery acquires a new- 
found significance. 

3.2 Elements of Photo Interpretation for Terrain Evaluation 

As it is not proper to call an imagery interpretation as 
photo interpretation, a new term visual interpretation of satellite 
imagery (VISI ) is coined for subsequent use. 

The visual interpretation of imagery for terrain evaluation 
is based on the spectral response of the ground objects* The key 
elements for the evaluation and interpretation are listed below. 
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a) Topography or ]a ndf orm 

j- 

b) Surface drainage pattern 

c) Spectral responses or gray levels 

d) Vegetation and Landuse 

All the above elements are studied sequentially using 

imageries of ! Bands 1, 2 and 4 of Landsat-4. 

3,2,1 Land. form 

Landform, in its broadest sense, implies the shape of the 
b nd i.e, topography* Each landform and bed rock type has its 
own characterstic topographic form, including a typical size and 
shape. In fact, there is of tec a distinct topographic change at 
the boundary between two different landforms, It is well known 
that once the origin rocktype- landform of a given area has been 
established, a large proportion of the interpretation is complete 
and the remainder will follow with comparative ease. The first 
part of the composite noun i.e. origin includes igneous, sedimentary 
and met am orphic and other fluvial origins. 

Under each of the primary origins applicable to rocks 
(igneous, sedimentary and metamorphic ) there are numerous different 
materials , These materials differ in many aspects; texture, 
mineralogical composition, nature of original source, secondary 
origin, etc. Based on these character sties, these materials may 
be classified into families or 'rock types'. 


The features such 
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as hills can be identified by their suggestive shapes. 

3,2.2 Drainage Pattern 

The drainage pattern and texture seen on any imagery are 
indicators of landform and bedrocktype and also suggest soil 
characterstics and site drainage conditions. Six of the most 
common drainage patterns are illustrated in Figure 3.1. 

The dendritic drainage pattern is a well integrated pattern 
formed by main stream with its tributaries branching and rebranching 
freely in all directions and occurs on relatively homogeneous 
materials. The rectangular drainage pattern is basically a dendr- 
itic pattern modified by structural bed rock control such that the 
tributaries meet at right angles. The trellis pattern consists of 
streams having one dominant direction and tributaries at right 
angles. The radial pattern is formed by streams that radiate 
outward from a central area as is typical of volcanoes and domes. 

The centripetal pattern is the reverse of radial drainage pattern. 
The deranged drainage pattern is a disordered pattern of aimlessly 
directed short streams, ponds, and wetland areas typical of 
ablation gracial till areas. 

Based on drainage texture, they are broadly divided as 
coarse textured and fine textured drainage patterns* The former 
develop where the soils and rocks have good internal drainage 
with little surface run off. It is reverse with the latter 
pattern and observed on easily eroded rocks such as shale. 




• f 

Figure 3 *4-' Six basic drainage patterns. 
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3.2.3 Spectral Responses 

The spectral response refers to the 'tone' or * brightness 1 
at any point on the imagery* As the satellites have sensors with 
wide ranging wavelength s, interpretation of the surface features 
becomes very easy. For example, the water bodies v;ill be darkest 
in tone (around 8-15 in grey level )in Band —4 imagery. But they 
can be better mapped from Band*~2 imagery where they appear somewhat 
lighter in tone against a darker back ground. Thus, varying 
spectral responses for various wave lengths can be made use of in 
visual interpretation. False coloured composites and Thermal JR 
imagery may help identify some more features. The interpreter 
has to develop his own keys while dealing with these wide 
ranging products, 

3.2.4 Vegetation and Land Use 

For identifying various vegetation types and changes the 
colour-lR film is proved to be the best one. With MSS imagery 
this process becomes almost impossible. Only vegetation cover 
we can detect is on hill tops. This identification is easy 
because of the peculier dark tone of the cover amidst its lighter 
surroundings. While searching for this element of interpretation. 
One should be knowledgeable about its changes through the year 
and accordingly he should evolve keys. IXaring the crop period, 
the true tonal patterns of soils are ‘hidden’ as only the crop 
is responsive to spectral analysis. 
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3*3 Ide nt If jcat ion of Roc k Types 

An effort is made to identify the various rock types 
pre sent in the area under study. After thorough VXSI of the 
three band, imageries (Bands 1, 2, and. 4 of Landsat *-4), it is 
concluded that the area primarily consists of two types of rocks 
viz, sedimentary and igneous * Again in sedimentary type, the 
consolidated and unconsolidated regions are clearly demarcated 
using B— 2 imagery. The boundaries . super at ing these rock types 
arc drawn v/ith exact precision. But, the bedding , jointing etc. 
ore not prominent as the scale of the imagery is too large. 

We first consider the characterstics of that helped 
identify the sandstone and other types latter, 

3.3,1 S andstone 

Topography ; massive, seemingly flat topped 

Drainage ? no drainage. The area is very small and seemingly 

un inhabited. No gullies because of excellent interval 
dr ainage 

Spectral response ; Light toned because of high quartz content. 

As it is surrounded by darker Recent Alluvium 
it {tan be mapped from both Band-1 and Band— 2 
imagery. 

Vegetation and Land Use: Sparse vegetation. Land use is rare. 

To confirm the above decision, a field test will do but 
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is not ncccssiatecl as the Geological Survey of India maps furnish 
all such details. 

3*3.2 Recent Alluvium (Deltas) 

Deltas form where streams discharge- into bodies of quic-t 
water. The delta area surrounding Krishna and Godavary ateas is 
distinctively seen on Band-2 imagery. This region is extremely 
flat interrupted only by irregularities associated with distribu- 
taries. Surf ace -drainage patterns in this area arc constructional 
i.o. the drainage patterns whose development has been modified in 
a pronounced, way by the operation of depositions!— geologic 
processes. Because of the huge sediments Carrie' by these rivers,, 
the region has become very fertile and most productive. 

3 *3.2.1 VISI of Delta s 

Topography : Nearly level surface bounded by upland areas and 
Water Krishna, delta has modified arcuate shape. 

Drainage and Erosion: Distributary streams present. The delta 
has manyrnajor channels arranged in a fan-shapec* 
pattern* Krishna River and other distributory channels 
are typically braided. 

Photo Tone: The soil is mdist because of the shallow water depths 
(0.5 m to 5 m). Hence low reflectance values (20 to 
30 in Band-2) are observed in the region* 
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Vegetation ancl Landuse : Extensively used for agriculture, ihe 

rock type in this region is identified 
as Recent Alluvium (Sedimentary - uncon- 
solidated) of fluvial origin. 

3.3.3 Salinity Survey 

The areas where high salt concentrations (saline or alkaline) 
are present, can be mapped with the help of MSS imagery, llic 
surface expression such as white encrustation will help interpret 
a part of saline area. In these places, crop vegetation may be 
extremely sparse or non-existant or may be composed largely of 
halpphytic plants. 

Determination of chemical composition ancl pH values of 
salt- impregnated areas is the work of field tests. 

White encrustations is clearly seen on Band-2 imagery 
and to some extent on Bancl«*4 imagery. The increase or decrease 
in salinity may be measured from imageries obtained periodically, 

3.3.4 Study of River Geometries 

The courses of major streams are displayed in unparalleled 
clarity and detail on satellite imagery of Band-2. Every curve, 
channel, clistributory is apparent. So, to study the horizontal 
geometries of a river, it is sjmply obtaining the lanclscencs in 
which the river of desired course is sensed and arranging them as 


mosaic 
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A frequent problem in unmapped areas, usually remote or 
relatively inaccessible, is the determination of an approximate 
stream profile. Satellite pictures provide, in such cases., a 
means d£ analysis to study the horizontal geometries of the stream 
under steady, 

3,3.5 Study of Surface Drainage 

In many cases, the most carefully prepared and large-scale 
topographic maps seldom show the tertiary and rill and gully 
drainage ways. Yet, these very unmapped drainage ways often cause 
much of the trouble for small development projects,. Also, it is 
those drainage channels that actually define the extent of any 
watershed, large or small. Consequently, their omission may be 
considered a serious flaw in many maps. 

Mapping the surface drainage net in all its detail from 
satellite imageries is easy, quick, relaible and economical. 
Against the darker background of sedimentary unconsolidated rocks, 
the water bodies, appearance is distinct. Either paper prints or 
a black and white film positive may be used for tracing than, 
Nagaraju (1986) rode on attempt to plot them. But qs that map 
does not contain all details like tertiarios, it is revised and 
presented as Figure 3.2, 
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3 # 6 I do nt if ic at 1 o n of U nc lassif iecl Rocks in the Study Ar ea 

The B aneW 2 imagery of the study area, as observed earlier, 
is excellent for soil as well as surface feature interpretation* 
Regarding the identification of sedimentary rocks (consolidated 
and unconsolidated), no serious effort may be made. But the rest 
of the area, identified as unclassified crystalline rocks by G,S,I*. 
presents a complex picture* It varies in tone unlike alluvial 
zone. Tliis region consists hills, ponds and other features wnieh 
are not there in the alluvial zone. This zone, as per GSI, 
consists both Igneous and Metamoridilc rocks which appears in 
light tone. 

Charnockite, that is. present adjacent to alluvial zone, 
is not distinguishable from the available three imageries. Tn.is 
may be due to its some common composition with its ncighbournood. 

This zone is shown white in the Figure 3*3 


3.7 Mapping of Built-Up Area 

For mapping of the built-up area in deltaic region, first, 
one- or two prominent towns are located on the image with the help 
of topo sheets . Iheir tone is grey and that of surroundings is 
dark. The towns and villages in this deltaic region appear as 
patches to specks depending on the extent of thoir coverage. This 
grey tone of built-up areas in deltaic zone is somewhat uniformly 


maintained 
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Mapping of this feature from unclassified crystalline 
rocks (mainly Gneisses) is extremely difficult. They are no 
longer conspicuous by their tone as the surroundings also more or 
less maintain the same tone. So, mapping from this rocktype is 
left out. Figure 3.4 shows the map of built up areas with some 
prominent towns and cities identified. 

3 . 8 Li near s 

Any linear feature in the landscape which possesses an 
abnormal degree of regularity is believed to be the surface 
expression of some structural feature in the bedrock. 

Linear s may represent one of several structural features; 
faulting/ bedding/ jointing, schistosity, gneissosity, contacts 
or narrow dikes. To identify a linear as one of these features, 
additional information from other sources is required. Hence, 
they remain as unidentified linears. 

The amount of detail which is desirable in plotting linears 
is a matter of experience and judgement. The chief trends and 
character of the linear pattern as shewn on ' 1 £ AGR^M 1 1 
(Fig. 3.6) must be shewn as well as its relative density. In tne 
initial phase of interpretation, or for a beginner, rt is probably 
better to err on the side of too much detail; the pattern can 
always be thinned out at the final interpretation if it tends to 

. The map showing the linears is 


obscure other information 
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labelled Figure 3*5, The linear s shown as broken lines along the 
shore line are may be due to sea transgression. 


3,9 Observations and Conclusions 

To sum up the salient points noted during the visual 
interpretation of satellite MSS imagery - 

a) Against the backdrop of alluvial soil, mapping of built up 
areas, both rural and urban/ is very effective ancl exact from 
Band-2 image, 

b) Sediment dispersion at river mouths can also be studied 

from Band-2 but Band-1 picture presents it in a better w ay. Band-4 
picture does not show any sediment discharge at river mouths 
because of the complete absorption (grey level falls down to a 
bare minimum of 4 ) of the near — XR electro magnetic energy* 

c) Locating the marshy land around the Kolleru Lake ancl its 

correct boundaries is not possible from Band-2 imagery, tne 
reason being both the marshy lane! and the surrounding alluvial 
soil appear dark in it. For this purpose, Band-4 picture is used. 

d) Locating small villages situated near the coast is difficult 
rather it is not possible, as they arc situated amidst forests 


and crosCkg 
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g) Tho hillocks can be traced better from Band-4 picture than 
those of Band 1 and 2, In band-4, their tone is light and are- 
conspicuous by their shape and size, 

f ) Surface water features are ma pped with remarkable ease in 

deltaic region fran Band-2 imagery. 
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CHAPTER 4 

FUTURE RECOMMENDATIONS 


With the whole range of satellite products, it is certain 
that VISI becomes very effective. Other features like vegetation 
that are not traced in this work will definitely be identifiable 
in false color composites where two or more bands are mixed to 
bring out some hidden features* 

Band ■- 1 picture of Landsat ^4 may be used to map the 
vegetation with the help of A d.ditive Color Viewer (ADDCOL), This 
requires 70 mm square positive films and the system* It is not 
carried out as both arc* not available with the Institute* 

The products of high resolution satellites may yield better 
results. For example, the SPOT simulator data contain a greater 
range of grey -'level values for all water areas than do the Landsat 
MSS data. The greater spartial resolution of the SPOT simulator 
data provides informations about small-scale hydrodynamics not 
available on Landsat MSS data, SPOT simulator data are found to 
delineate water masses with a high degree of seperation 
(STEVEN ot al, 1985), The SPOT data are proved much useful for 
agriculture, evaluation of geological alterations, forest cover. 

With SPOT or Thematic mapper data, the K -Class Classifier 
results may be highly reliable because of high degree of resciluticm 
It may be interesting to watch how the classifier reacts with] two 
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banc's as input signals, 

A soil, c 1 assi £ ic at io n map showing all the three classes of 
the study area may be prepared with the help of classifier, The 
area near Guntur where all the three present may be chosen for the 
purpose. Similar ly, Using the classifier , urban areas and their 
growth can be mapped. Another important application may be flood 
mapping. In fact, for flood mapping and show-emver mapping and 
a simpler method called Density Slicing may be employed. 

As the ground water tabl^ depths of the study area are 
available with the Department, it can be syfh with the help of 
classifier how many classes are present and also a computer map 
showing those divisions can be prepared, 

Ihe results of this classifier may be compared with those 
of other classifiers such as Bayes’ Classifier, 
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:omrf»rinds tor a tunnel CHANGE program 


A SR MTA;* (,;(, 
r change 
M A K t'l 


HKl’ftJ »> 

O.SK : ¥ . *-HTA« ; t q t/UKc : i J a3/oi;tf : t oOoahjK : 


10/MODk: EBCDIC/ NU ERROR/NOCKLiF 


1 f: h P C H A N c; E T [\ X T 


lommands : 


) AT A 

lAKf'J 
:xix 
IELP 

i n: 

ttCTATN 
ion 
:rask 

'RTNT 
twitches 


the !u\,iie of the data file to he used as in 
a normal command string ,, MTA«nuV:FiLE.EXHt , *PNJ , 
CaUTAUl.K.hAV to create tine conversion tables. 
Terminate the program. 

This little message, 

Ca-U 1.1 IK, HUT into core and execute It. 

Allows the user to accumulate commands, 
perform the current command string. 

Erase retained commands. 

Print, tne current command string. 


!ufters;x Use x buffers. 

idvance:x Advance x files before operation, 

lacK.space;x Backspace the tape x files before operation. 

Slock: x set blocking factor to x, 

iecord;x Set record size or size of largest block to x. 

>ensity:ard Set mag-tape density to am targ=zoo ,55b , BOOT „ 

letain Retain the following commands, 

lun Perform tne current command and retain it, 

lelp These, few hints. ^ s . , 

labelsftry Set label type as described below for mas-tape, 

iodesary Set file character set as described below, 

>arity:arq Ret maq-tape parity [oda-odd. even^even, def aul t=odcil , 

•asswora ; arg Set tue password for GE .labels. . 

ieel:x set serial number in labels that nave this feature,, 

£ industry Initialise for industry compatible 9-channel tape, 

: Scan scan tape for file named. 

5 Error Don't ignore checksum and parity errors. /(11 

< Span Records cross blocks [Implied if "/block: O' ] . 

c Rewindsarq Rowing tapes Card-before, after, always, omitcJ 
‘ Unload r Unload tape after operation. 

< Tell Type tile names on a wild card search. 

‘ List last the device directory, J . _ 

! FUst last the device directory [file names only!. 
c Header Print headers on the line printer, 

1 Crlf Ascii file nas crif's, 

lots: To turn a sdtcn. off concatenate "no” with the switch, 

Switches flawed with an (*) have this feature. 

'babel - switch modifiers; 


one Maq-tape has no labels just data [default], 
ine Special labels for pECfiyatew-lO. Tah ^ 1<= 
iqital Process labels as standard DIGITAL labflf . . _ file 

urroudhs Process labels as standard burroughs lapels, 

BM Process labels as standard JgW labels . 

E 6 3 5 Process labels as standard GE-635 labels. 


Change Text contd, . 


Jot.e: fifths .1 s are written in the mode specified except for IBM 

nu 1 GE-bJB lapels. IBP lapels are written in BCD lor 
'-track drives and KDCIHC tor 9- track drives, Cr J * , b35 
labels are always written in Gf>BCn. 

'tfode « switch modifiers* 

\SCLl Pile character set: is 7-hit ASCII. 

U'ASCXI Pi 1 e. character set is H-bit ASCII, 

JiSAiSCII Pile character set is 9 -bit ASCII, 

in AGP Pile is read and written as 36-bit words* 

> 1 X B I T File character set is SIXBTT, 

r I X s T. x Pile character .set is SIX BIT with no control words, 

iCh File character set is BCL. 

)C7> Pile character set. xs BCD. 

5KBCD File character set; is GS>BCL), 

10NBCD File character set is HON^YwELli BCD* 

SBCDIC File character set is fixed FBCDIC, 

MSBCIUC File character set is variable tlBCDlC, 


Phe end of the input record is de termined by the mode of: the 
Input file as described below: 


\SCII 

5IXBTT 

axsix 


\C\, 


1CB 


BCD TC 


ecord 
BCD IC 


ecord. 


Terminated by a carriage-return character or record count, 
Peter mln ted by the header word in front of each record. 

Always copies tne number of bytes specified by the. record size 

Always copies the number of bytes specified by the. record size 

Always copies tne number of bytes specified by the record size 

rpixedJ Always copies the number of bytes specified by the 

i £* m 

[Variable]! Determined by the header word in front of each 



Cn general all commands may be abbreviated to any number 
3 f characters that; win allow that command to be unique. However 
io more then six characters are checked for validity in any of 
the commands. 


rhe characters "?" and "* H may be. used to denote a wild card tor files 
on mag- tape f disk f or dec tape. The character "P" denotes a 
command file which change will read for commands* If only 
an aitmode is typed change will enter dialog mode. 



APPENDIX - 2 


Sof tware 

Under this head, the various programs used for K— Class 
classification and mapping are explained. The listings arc 
reproduced in the latter pages. 

1 , K— Glass .FQ _R 

This program is for the K --Class classification algorithm 
developed by Zagalsky (1968) which gives the results in the form 
of a confusion matrix. This program developed by Serreyn and 
Nelson (1973) is modified slightly to apply for the present study. 

Explanati on of Various FORTRAN Variables of the Program 
NCLASS ~ Number of classes of data 

NSIG - Number of signals supplied for each feature 

NFEAT - Number of features or Bands 

P (I) ~ F is the apriori variability of occurence and P(l) is 

the probability of occurrence of class I. 

B (I, J) ~ B is a matrix multiplier of the feature vector X* It is 
calculated in the training of the classifier. 

C (I) - C is a vector of constants calculated in the training 

of the classifier. 

Cov (1,1) ~ Cov is the covariance matrix of attributes 

AINCOV (1,1) - Inverse of the covariance matrix Cov. Inverse is 

calcualted by calling the subroutine MATINV . 



E (i) - E is an interim vector variable 

D (I) “ D is thG decision vector and D (I) is its Ith clement. 

NOCL (I) - is the cumulative distribution of the data NOCI, (2) 
is the number of data for class one + class two.. 
Label (i) - it is to label tte output of the classifier 


i . OB — It generally equals 1 since new data of different features 
arc to be used. NPROB # 2 is used if we have unequal 
samples for each class but wish to have equal a priori 
probabilities for each class. 

IDEC (i/J) — IDEC is the confusion matrix with the correct 

identifications as the diogonal elements and the 
incorrect ones as the other elements of the matrix. 


It is valuable to note that the B and C matrices are not 
affected by the a priori variabilities of the classes; they arc 
based upon the sample data only. Another valuable item to note 
in general, which can be used to check the computer program is 
th at 


N class 

B (I,J) P(I) = 0 

1=1 
N class 

c (i) p (i) = o 

i=i 


These equations are baraed upon the theory that sum of the D (l)'s 
must equal one since the sum of P (X)'s equal one^ 



Inpu t 


'Ihfc input for this program is supplied from two files; 

FOR 24 .DAT and F0R48.DAT. FOR24.DAT consists the following four 
data cards, 

10 Label (I) , i = 1,20 
20 Labcll (I), I = 1/20 

30 Nclass / Nfeat Nsiq. 

40 NOCL (I), I = l / Nclass 

FOR 48, DAT contains unformatted video data. The number 
of data should coincide with NSIG of FOR 2 4, DAT, 

O utpu t 

The output file is FOR 52. DAT, Besides the formatted 
input data/ it consists of decision table, confusion matrix, 
percentage matrix for the classification and decision boundaries. 
Boundaries among the various classes is calculated by subroutine 
BOUND in K-dimo nsional space where K is number of classes, 

2. KCIMAP ■ FOR 

This program is for generation of a computer map based on 
K-Class classifier decisions. The area to be classified may be 
located from TOPO sheets. The co-ordinates of this area arc to be 
converted into scan line numbers and pixel numbers. To have an 
idea about the number of classes present in this area may be had 



from topo sheets or paper print imagery* The features suitable 
for the classification may be decided by feeding in the sample 
data in the program KCLAS s. FOR. The combination oiff bn ride ? *thrrb .gives 
maximum correct identif ication is found * After deciding the 
number of classes present in the area, and the suitable bands for 
classif ication, the data from the entire area to be classified is 
retrieved from CCT and stored in individual files. The program 
KCLMAP-FOR allows two bands and .9»,6oo signals per each bands. 

The K-Class classifier decisions for individual signals 
are assigned K-diffcrent symbols for easy identification on the 
map. These symbols are to be discreetly chosen so that the 
various classes are distincly visible. Those symbols are 
supplied to the program through an array GREY declared as an 
integer subscripted variable. After assigning these symbols to 
the classifier decisions, they are stored in the array JAR. If 
it is desired to use higher number of signals, the capacity of 
JAR may suitably be altered. 

It is to be remembered that the width of the physical page 
is 132 characters. Hence, if the width of the area to be 
classified i.e. difference between end pixel and starting pixel 
exceeds 132, the program demands amendment. By dividing 
the area into the two parts and classifying individually wc can 
by pass this amendment. But, Classification of the area as one 

unit givos beist results. 



3 . Record .FCR 


This program is to read a scan line of a scene. Each 
line is a record on the CCT. CCT contains 3 files viz. Tape 
Directory, Teipe Header and Video Data, If nth line is to be 
retrieved, the first two files and 4(n— 1) records one to oe 
skipped. This is because the tape is in BIL (Band Inter l/oavec. ) 
form at . 

This program reads a record of CCT and then processes tno 
data (from Binary to ASCII ) . It picks up 36 bits at a time 
(a WORD) and then it divides into 4 bytes (pixels). Ihe divisive 
constants are as follows t 


bits/bytes: 0 

constants :: 



2**4 


2**12 


2**20 


35 

2**28 


This processed data are reflectance levels that range 
from 0 to 255. If a particular grey level is shown negative, the 
true level is equal to(256 + grey level). 

The divisive constants discussed above are valid for 36-bit 

„ , . n1?r eW j flC! por 32*-bit word machines, the 

word machines like DEC series. 

,, , „ 2**16/ 2**8 and 2**0 

constants will' correspondingly b / 

(i*e. 1 ). 



4* PIXEL *FOR 


If a good number of grey levels are to be noted in a 
record then the above program may be used. On the other hand, if 
the level of one pixel is to be noted, the program PIXEL* FOR may 
be used. Ibis is a modified program of RECORD.FOR. After 
execution of this program, the user has to type in the pixel 
number whose grey level is required. The output is from the file 
PIXEL. DAI. The output also contains grey levels of five proceeding 
and five succeeding pixels of the current ono. This program 
counts the length of initial zero fill and adds it to the pixel 
number supplied, 

5 * PART. FOR 

For classification purpose, one may have to retrieve a 
part of the scene*. Tor this purpose, this program may be used 
to retrieve a part of the record. A file with extension MIC 
may be used for storing the data and execution of the program 
if it is desired to retrieve a large number of records. 

6. DEN&TY . FOR 

This FORTRAN program is for Density slicing, a simple 
.method for classifying the data based on their grey levels. If 
the range of the' grey levels for a particular class is definitely 
known, then this method is very effective though it is a primitive 


method 



'Ihis Density Slicing method is synonymous to HIGHLIGHTING 
in Image Processing where a contour level and Interval is fed.- in 
to see the highlighted feature on the video screen, Ihc program 
DENSTY.POR is an interactive program where under execution the 
user has to supply the number of classes he intends to classify/ 
the characters to represent them, contour level and interval lor 
each class and size of the input matrix sequentially. The output 
is stored in Map 1, Dat to Map 8, Dat , As cautioned previously, 
the user has to bear in mind that the maximum width of the output 
from line printer is 132. So, if the width of input matrix exceeds 
this number , more than one output files are to be specif led. 
According to this program, if the width of the matrix is 480 f the 
user has to specify that 480/120 = 4 files are required for output. 
Tliis is done during interaction, 

7, LINPIX, FOR 

The ground co-ordinates, latitude and longitude of a point 
will be substituted by scan Zinc number and pixel number on the 
imagery. So, for finding out the corresponding line and pixel 
number of a particular ground point, its latitude and longitude 
ore to be precisely noted from topographic maps. These are fed 
in as input to get the scan line and pixel number. The program 
consists constants obtained ffom solving the eight equations that 
have six unknowns. These const arts are valid only for a particular 


scene 



8. SORT. PAS 


'Ihe line number anc! pixel numbers obtained from above 
program are to be processed before ret riving the data i.o. they 
are to bo arranged in ascending order in order to facilitate the 
retroival. This program also gives the number of records to be 
skipped at each stage, the scan lino number division and Reaord 
number division. This data will be extremely helpful when 
particularly a large number of records are to be retrieved. The 
scan line division, Record number division obtained from the 
program should tally with those of the retrieved one. 

9, PLOT, FOR 

'Ihis is a graphics program written in FORTRAN that uses 
General Purpose Graphic System (GPGS) subroutines, for generating 
graphs or scatter diagrams. 

1°. FIG, PAS 

This is a pascal program for generation of a line printer 
picture of an area by feeding in the grey levels as input,. The 
output will be in 32 levels whereas the input data will be in 
256 levels of intensity. The characters for representing the .# 
levels is a difficult task and so nc-cds careful design. 
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